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Abstract 



Making use of the effective field theory expansion recently developed by the 
authors, we compute the electromagnetic form factors of the deuteron ana- 
lytically to next-to-leading order (NLO). The computation is rather simple, 
and involves calculating several Feynman diagrams, using dimensional regu- 
larization. The results agree well with data and indicate that the expansion 
is converging. They do not suffer from any ambiguities arising from off-shell 
versus on-shell amplitudes. 
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I. INTRODUCTION 



The techniques introduced by the authors in refs. pj put the study of low energy two- 
nucleon interactions on the same footing as chiral perturbation theory in the mesonic and 
single nucleon sectors @. In particular, there is a systematic low momentum expansion, such 
that at any given order one need only calculate a finite number of Feynman diagrams to arrive 
at an analytic result. The procedure is superior in several ways to the conventional technique 
of solving the Schrodinger equation with a potential constructed to fit the scattering data: 
(i) There is a well defined expansion parameter, and one can estimate errors at any given 
order in the expansion; (ii) it is straightforward to incorporate relativistic and inelastic 
effects within the expansion; (iii) analytic results allow one to see quite simply the relative 
importance of short- and long-distance physics to a given process; (iv) there is no ambiguity 
concerning off-shell matrix elements when calculating physical processes; (v) at low orders 
in the expansion, the number of free parameters to be fit to the data is few, and the same 
parameters are used in all processes. The results at lower orders in the expansion are 
therefore very constrained. 

Until now, the techniques of ref. [|IJ have only been applied to reproducing scattering 
phase shifts. While a necessary first step, fitting the phase shifts does not seriously test 
the method, as the low energy phase shifts can be well fit by rather simple functions of few 
parameters. What is needed are calculations of dynamical processes that involve the same 
interactions as are fit to the iViV phase shifts. The obvious ones to consider are NN — > NN'y, 
np — > c/7, parity and isospin violation in NN processes, pp — > de + u, and the deuteron 
electromagnetic form factors. In this paper we present the perturbative calculation of the 
deuteron electromagnetic form factors at NLO. This subject has been addressed previously 
in the context of effective field theory in refs. [[|f|], although using a somewhat different 
formalism and involved numerical, as opposed to analytical, calculations. We preface the 
calculation with a brief review of our expansion, and a discussion of the deuteron. After 
identifying the graphs contributing the electromagnetic form factors, we show explicitly that 
there is no ambiguity arising from the fact that the nucleons in a deuteron are not on their 
mass shell, even though the couplings in the effective theory are fit to iViV scattering data. 
We conclude with a discussion of features that will appear in the NNLO (next-to-next-to 
leading order) calculation of the form factors. 

It is not unreasonable to ask why it is worth pursuing an effective field theory descrip- 
tion of the deuteron since effective range theory can be used to predict many of its 
properties |7|||. For some quantities, like the deuteron charge radius, effective range theory 
is remarkably precise. The primary motivation is to make a clear connection to QCD and 
therefore enable systematic calculations to be performed, even for processes where effective 
range theory is not applicable. Furthermore we expect that, even in cases where effective 
range theory is very accurate, the effective field theory approach will surpass this level of 
precision if pursued to higher orders. 

In this paper the application of the effective field theory expansion of ref. || to processes 
involving the deuteron is developed. For definiteness we focus on the electromagnetic form 
factors of the deuteron. However, it is straightforward to use the methods developed here 
for other quantities, like the cross sections for np — > and 7c? — > 'yd. 

At NLO the predictions of the effective field theory expansion for the electromagnetic 
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form factors of the deuteron are not as accurate as those of effective range theory. However 
at NNLO they should reach the precision of effective range theory. Furthermore, the effective 
field theory approach is a systematic one at some order it will include physical effects beyond 
those that are encorporated into effective range theory. 



II. EFFECTIVE FIELD THEORY FOR NN INTERACTIONS 

In order to compute the electromagnetic form factors of the deuteron, we must consider 
the possible interactions between nucleons, pions and photons. In an effective field theory, 
these interactions take the form of local operators, constrained only by the symmetries of 
QCD and QED. In this section we discuss the form of the operators that occur to the 
order that we will be working, and then turn to the issue of power counting, which allows a 
consistent expansion of the form factors. 



A. Interactions 

Terms in the effective Lagrangian describing the interactions between nucleons, pions 
and photons can be classified by the number of nucleon fields that appear in them. It is 
convenient to write 

C = C + A + C2 + (2.1) 

where C n contains n-body nucleon operators. 

C is constructed from the photon field N 1 = (A , A) and the pion fields II; it does not 
contain any nucleon fields. The pion fields are incorporated in a special unitary matrix, 

S = exp— , n=^_ ^ j , (2.2) 

where / = 132 MeV is the pion decay constant. E transforms under the global SU(2) L x 
SU(2) R and U(l) ern gauge symmetries as 

£ -> LTiR\ £ -> e ^-Se^ aQem (2.3) 

where L e SU(2) L , R G SU(2) R and Q em is the charge matrix, 

9- = (JS). (2.4) 

The part of the Lagrange density with no nucleon fields is 

C = i(E 2 - B 2 ) + ^Tr DpEDftf + ^cuTr m,(E + E f ) + . . . . (2.5) 
2 o 4 

The ellipsis denotes operators with more covariant derivatives D^, insertions of the quark 
mass matrix, m q = diag(m u ,md), or factors of the electric and magnetic fields. Acting on 
E the covariant derivative is 
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= + ie[Q em , E] A M , (2.6) 

The parameter u has dimensions of mass and m% = uj(m u + rrid). 

When describing pion-nucleon interactions it is convenient to introduce the field £ = 
expill// = Under SU(2) L x SU(2) R transformations 

£ ^ L£[/t = £7^, (2.7) 

where [/ is a complicated nonlinear function of L, R and the pion fields themselves. Since U 
depends on the pion fields it has spacetime dependence. The nucleon fields are introduced 
as a doublet of spin 1/2 fields 

N=(l), (2.8) 

that transforms under chiral SU(2)l x SU(2)r symmetry as N —>■ UN and under U(l) 
gauge transformations as N —>■ e tot< ^ em N. Acting on nucleon fields the covariant derivative is 

D^N = (fy + V^ + leQemA^N, (2.9) 

where 

V l > = \{tD l £ + $D ll £) . (2.10) 

The covariant derivative of N transforms in the same way as N under SU(2)i x SU(2)r 
transformations (i.e., D^N — > UD^N) and under U(l) gauge transformations (i.e., D^N — > 

The one-body terms in the Lagrange density are 

£ x = iV f (*A> + ^) N + '-f^a • (£D£t - ^D^A 



^ f «o + -^ICt 3 Z + ^r 3 ^] U • BN + . . . . (2.11) 



2M V 2 

To the order to which we are working k = §(/%> + K n ) and isoscalar and 

isovector nucleon magnetic moments in nuclear magnetons, with 

k p = 2.79285 , K n = -1.91304 , (2.12) 

at tree-level. At higher orders there will be contributions to eq. ( 2.12|) from pion loop 
graphs ||. The isoscalar magnetic moment k,q receives leading corrections of the form 
m^.log(m^/A£), suppressed by two powers of the pion mass. In contrast, the isovector 
magnetic moment k% receives leading corrections of the form m^, suppressed by only one 
power of the pion mass. The ellipsis in eq. ( |2.11 ) denotes higher order terms that do not 
contribute at the order we are working. 

Finally it remains to consider the two body operators. Some of these were discussed in 
refs. [[IJ ; however, since we will be computing electromagnetic form factors of the deuteron 
there are additional considerations that didn't arise in the NLO calculation of nucleon phase 
shifts. 
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First we will consider the two-body operators involving nucleons alone, then we will look 
at those containing a photon; to the order we will be working, we need not consider two- 
body operators involving pion fields. In the spin triplet channel, there is one NN contact 
interaction with no derivatives or insertions of the quark mass matrix, corresponding to a 
diagonal transition 3 Si — > 3 «Si; the coefficient of this operator is taken to be Co. There is 
an additional contact interaction involving no derivatives and one insertion of the quark 
mass matrix, with coefficient D 2 ; it can be distinguished from the Co interaction by its 
chiral properties. There are five contact interactions involving two gradients, corresponding 
to diagonal transitions in the 3 Si, 3 Po, 3 Pi, and 3 P2 partial waves, as well as an off-diagonal 
3 S\ — 3 Di transition. Only the first and the last of these are relevant for the deuteron; 
furthermore, at NLO we can ignore the 3 S\ — 3 D\ transition interaction. Thus the only V 2 
two-body contact interaction we will consider is 3 5i — > 3 S\, and has coupling Ci- Therefore 
for a 3 S\ — > 3 S\ scattering process, where the incoming nucleons have momenta pi, P2 and 
polarization i, and scatter into states with momenta p^, p 2 and polarization j, the Born 
amplitude arising from the contact interactions is 



iA 



-iSij 



C + D 2 ml 



C 2 /, 



p 2 r + 



P2) 



(2.13) 



The form of the C 2 amplitude is fixed by Lorentz invariance (which is equivalent to Galilean 
invariance to the order we work), and by the normalization we used in ref. [0, where in the 



center of mass frame, where we defined the amplitude to be —iC 2 p, p = |pi 



B As 



discussed in appendix [B], while one can construct a two-body contact interaction with one 
factor of do instead of two gradients, for any ^-matrix element (including those involving the 
deuteron) one can use the equations of motion to eliminate time derivatives for gradients. 
Thus no independent do contact interaction needs to be introduced. 

Including gauge fields introduces several two-body contributions to the electromagnetic 
current. Firstly, the C 2 interaction described above becomes gauged. Secondly there are 
two new two-body magnetic moment type interactions. In order to write C 2 compactly we 
define the matrix P, which projects onto the 3 S± state, 



Pi = —?=o- 2 aiT 2 



TrPM, 



-^p}pj = h 



(2.14) 



where the a matrices act on the nucleon spin indices, while the r matrices act on isospin 
indices. Then the two-body Lagrangian may be written as 



C 2 = - (C + D^Ttm,) (N T P,N)\N T PiN) 



(N T PiNy (n t pr5 2 + £) 2 p - 2r5pE> 



N 



+eL 2 \(N T PiN)\N T P i( r ■ BiV) + h.c. 



+ ... 



h.c. 



(2.15) 



1 The couplings Co, D 2 and C 2 are the same couplings that appear as Cg Sl \ D2 Sl \ and C 2 Sl ^ in 
refs. |Q; we drop the 3 Si designation here as there can be no confusion with analogous couplings 
in the ^o channel. 
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where the ellipsis refers to contact interactions irrelevant for the deuteron channel, or of 
higher order than we will be considering. The new coupling L 2 corresponds to an interaction 
that did not enter the calculation of NN scattering, but which affects the deuteron magnetic 
form factor. As written, eq. fl2.15|) is not chirally invariant, which can be remedied by an 
appropriate insertion of the £ fields; however, since the two-body operators with pions do 
not contribute at NLO, we omit them. 



B. Power counting 

We begin by summarizing the results of refs. [1]]. The starting point is the effective 
Lagrangian for nucleons, pions and photons introduced in the previous section. The part 
of the Lagrangian describing purely mesonic interactions, as well as interactions between 
mesons and a single baryon, is the conventional chiral Lagrangian. In addition there are 
local interactions corresponding to short distance interactions between two nucleons. These 
contact interactions are expanded in powers of derivatives and insertions of the quark mass 
matrix, m q . (Isospin violation from the difference between the up and down quark masses 
is neglected. Consequently insertions of m q are equivalent to factors of rn^.) The lowest 
dimension operator is a four fermion contact interaction; there are two independent operators 
of this form, corresponding to the x Sq and 3 S\ channels. The next lowest dimension two- 
body operators involve a factor of p 2 , where p is the momentum of one of the nucleons 
in the center of mass frame, or a factor of m 2 . There are seven independent p 2 operators 
corresponding to diagonal matrix elements in the {^o, ^1, 3 Si, 3 Po, 3 ^i, ^2} channels, as well 
as a 3 Si — 3 Dx mixing term; there are two independent m 2 operators corresponding to the 
^So and 3 5i channels. At higher powers of derivatives, the number of contact interactions 
quickly grows. 

Central to effective field theory is a power counting scheme which allows one to calculate 
consistently to any given order in the low energy expansion. A main point in refs. Jl| was to 
develop the PDS subtraction scheme which allows one to readily identify the order of any 
particular Feynman graph. The scheme involves computing loop diagrams using dimensional 
regularization, and then subtracting off the poles in dimensions D < 4, which correspond to 
logarithmic or power-law divergences. A typical integral in this scheme is 

/„ = i(p/2)*- D J -^L q 2 " (E/2 + <,„ - q 2 /2A/ + fe)~' (B/2 - <,„ - q 2 /2M + is) 



\A-D 



-M(ME) n (—ME - isr-»»r Jp-_ 1)/2 



PDS 



M 



^ -{ME) n [-) (ja - V-ME - is) . (2.16) 

The last step includes the finite subtraction mandated in the PDS scheme. The parameter 
/1 is the renormalization scale and physical observables are independent of it. In fact, one 
may set /1 to zero and recover the usual minimal subtraction scheme (MS) with fi = if 
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one wishes f\ However, a change in \i must be compensated by the renormalization group 
flow of the couplings in the theory. Therefore, what is a weak coupling at one value of \i 
can be strong at another, which effects how one defines the power counting scheme. 

Rapid scaling with // is only an issue for two body operators, and then only for those 
affected by the large scattering lengths in the ^So and 3 Si channels. Consider a four nucleon 
contact interaction connecting angular momentum states L and L', where conservation of 
angular momentum and parity requires \L — L'\ to equal zero or two. We assume that the 
operator involves m insertions of the quark mass matrix, and 2d = (L + L' + 2n) spatial 
gradients, and has a coefficient C!^ L n . By examining the coupled renormalization group 
equations in the PDS scheme, one can determine that these couplings scale as 

Lm ' nW otherwise. 1 ] 

in the region 1/a <C A^tv, where 

AnN = ~ 300 MeV ' (2 ' 18) 

Here M is the nucleon mass, qa = 1-25 is the axial current coupling and / = 132 MeV is 
the pion decay constant. Thus in the deuteron channel, Co ~ while Ci and D 2 scale as 
fi~ 2 . Extending the analysis to include photons, we find L 2 ~ /i~ 2 as well. 

The coefficients of the four-nucleon contact terms that have explicit factors of the elec- 
tric field E or the magnetic field B scale similarly to those in eq. ( |2.17j ), counting gauge 
fields as derivatives. For example, the L 2 operator in eq. ( |2.15|) counts as a two-derivative, 



L = L' = operator, and its coefficient scales as L 2 ~ fi . The rapid scaling of the oper- 
ators contributing to 5*-wave processes is what makes our expansion different than the one 
proposed by Weinberg JTO . 



Armed with the above results, we are able to arrive at a particularly simple set of rules 
for determining the order of a graph. Choosing the scale /i ~ p ~ m n ~ Q we perform an 
expansion in Q, where 

1. Each nucleon or pion propagator scales as Q~ 2 ; 

2. Each loop integration / d A q scales as Q 5 ; 

3. A gradient at a vertex contributes Q 1 , while each time derivative scales as Q 2 ; 

4. An insertion of the quark mass matrix m q at a vertex counts as Q 2 ; 

5. The coefficient of the contact interactions scale according to eq. ( 2.17| ). 



2 In the MS scheme with /i = one must first integrate out the pion to avoid factors that diverge 
as log(m 2 /[i 2 ). 
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The first three rules follow simply from the scaling of four momenta appropriate to the 
nonrelativistic regime. Explicitly, Mq ~ q 2 ~ Q 2 . The fourth rule is familiar from conven- 
tional chiral perturbation theory, m q ~ m 2 ~ Q 2 . The procedure for calculating physical 
quantities of interest is to write down the most general effective field theory consistent with 
gauge invariance, chiral symmetry and Lorentz invariance [], and then compute the desired 
matrix element to a given order in the Q expansion, following the above rules. Note that ac- 
cording to the power counting rules, a loop with two propagators entails a factor of Q, while 
the coefficient of the lowest order NN contact interaction (Cf m with L — V — d — m — 0, 
defined to be C ) scales as l/Q; thus any graph may be dressed by an infinite bubble chain 
with Cq interactions without changing the order of the graph. 



III. THE DEUTERON FORM FACTORS 

A deuteron with four-momentum p^ and polarization vector & is described by the state 
|p, e), where the polarization vector satisfies p^ = 0. An orthonormal basis of polarization 
vectors ef satisfies 
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i=l 



TjT, 



(3.1) 



where is the deuteron mass. It is convenient to choose the basis polarization vectors 
so that in the deuteron rest frame ef = 5f . Deuteron states with these polarizations are 
denoted by |p, i) (i.e., |p, i) = |p, ef}) and satisfy the normalization condition 



(p',j|p^) = ^(2vr) 3 5 3 (p-p')^ 



(3.2) 



In terms of these states and to leading order in the nonrelativistic expansion, the matrix 
element of the electromagnetic current is 



(p',j\Jem\P^) = e 



2M, 



F c{q 2 )$ij + ^| F e(<? 2 ) (w - 

F c (q 2 )5 tJ (p + p') k + F M (q 2 ) - Sfa 

1 „ , / 1 



2MV ^ ' 3 



F Q (q 2 ) cuqj - -q 2 ^ (p + P' 



(3.3) 



where q = p' — p and q = |q|. These dimensionless form factors are normalized such that 



3 Relativistic corrections are accounted for as perturbations according to the above power counting 
rules, and at the order we work the theory only appears Galilean invariant. The procedure for 
dealing with relativistic corrections perturbatively requires distinguishing between potential and 
radiation pions at NNLO, as discussed in [|J. 



S 



FIG. 1. The perturbative expansion o/E. The first row has the leading 0(Q) result; (g) repre- 
sents an insertion of the interpolating field defined in eq. ( \3. % ). The second row has the complete 
subleading 0(Q 2 ) contribution, where ■ and o denote the C2 and D2 interactions respectively. The 
third row shows a couple of 0{Q ?J ) NNLO contributions, which we do not calculate here: the ex- 
change of two potential pions, and the dressing of Co (the pointlike NN vertex) by a radiation 
pion. 



2M d 
1 



F c (0) = 1 , 
Fm{0) = Hm 



m: 



(3.4) 



where /jm — 0.85741(e/2M) is the deuteron magnetic moment, and /ig = 0.2859 fm 2 is the 
deuteron quadrupole moment. 

As shown in appendix the form factors are readily calculated by computing in pertur- 
bation theory the irreducible two-point function £, and the irreducible three-point function 
r^ 4 . In the present context, "irreducible" means the sum of graphs which do not fall apart 
when cut at any Cq vertex. The matrix element of the electromagnetic current is then given 
by the exact relation 



\P',jVemlP^> 



£|(£Xq) 
cAY,(E)/cAE 



(3.5) 



E, E - 



where B is the deuteron binding energy and E is the energy of the incoming two nucleon 
state in the center of mass frame, 



E = E 



AM 



+ ... 



E = (p°-2M), 



(3.6) 



where the ellipsis refers to relativistic corrections to the energy-momentum relation. E is 
the analogous quantity for the outgoing nucleon pair. By Lorentz invariance, S and r M can 
only depend on the energy and momentum in this combination. 

We can now expand the relation eq. ( |3.5|) in perturbation theory and determine the form 
factors by comparing the result with eq. ( ^.3[ ). The two-point function has the graphical 
expansion shown in Fig. [1], where the ® vertices represent the insertion of an interpolating 
field T>i with the quantum numbers of a deuteron with polarization i. We take T>i to be 
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FIG. 2. The expansion of T°. In all of these graphs, the photon corresponds to Aq with the 
minimal coupling to the proton propagator, arising from the gauged nucleon kinetic energy term. 
The graph in the first row is the leading 0(Q~ 1 ) contribution, ^9^. The second row are the 
subleading graphs at O(Q ), summing to give In the third row are several graphs contributing 
at the 0(Q l ): a dressing of the photon-nucleon vertex, a relativistic correction to the nucleon 
propagator, and an exchange current contribution. 



V t = N T P t N, (3.7) 



where Pi is the projection defined in eq. ( 2.14j) . The form factor one calculates does not 



depend on the particular choice for T>i, so long as it is used consistently. 

By examining the graphs and using the power counting outlined in the previous section, 
one sees that S begins at order Q 1 — the leading graph has two nucleon propagators and 
one loop. At subleading order, 0(Q 2 ), there are three two-loop graphs, one involving the 
exchange of a potential pion (which has a derivative coupling), one with an insertion of 
the C 2 (D) 2 two-body operator, and one with an insertion of the D 2 m 2 two-body operator. 
Recall that with renormalization scale \i ~ Q the coefficients C 2 and D 2 are 0(Q~ 2 ). At 
0(Q 3 ) there are a host of diagrams, including the exchange of two potential pions, or one 
radiative pion, as well as p 4 relativistic corrections to the nucleon propagator, etc. We have 
calculated £ to 0(Q 2 ), and the results are presented in appendix [A[ 



A. The NLO computation of the electric form factors 



To compute the electric form factors F c and Fq we need to calculate the three-point 
function T^-, which is expanded graphically in Fig. 0. The results for both the leading 0(Q~ l ) 
and subleading 0(Q°) contributions are presented in appendix 0. Once T is computed in 
the Q-expansion, the electric form factors can be determined by expanding eq. {\3.i 



as 



(p',j I^LIp,*) 



F° 

M-i) 



d£ ( i)/d£ 



+ i 



r° 0) (d£ (1) /d£ 



F° 

1 (-1) 



d£ (2) /d£ 



\ 2 



d£ ( i)/d£ 



+ 0(Q 2 ), (3.8) 



where T^, £( n ) denote the 0(Q n ) contribution to T and £ respectively. We have suppressed 
the q dependence of T , and its polarization indices. Furthermore everything is evaluated 



on-shell, E = E = -B. Since d/d£ ~ 0(Q~ 2 ), the first bracket in eq. (p]) is 0(Q°), the 
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second bracket is 0{Q 1 ), etc. Therefore, taking into account the explicit factors of q in the 
definition of the form factors, eq. (|3.3| ), we see the electric form factors have a Q expansion 
of the form 

<o) , p(l) 



a +F£> + 0(Q< 



F Q = F { Q 2) + F { Q 1) + 0(Q", 

where ~ 0(Q n ). 

Using eqs. ( |3.8| , [7^6| , |Al7| ) gives our leading result for the electric form factors, 



(3.9) 



4V 
Q 



4 7 f -i 
— tan 

Q 

, 



47, 



where we have defined 
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(3.10) 



(3.11) 



The subleading form factors are extracted from eqs. ( |3.8| . |A6|fA18| ), and presented in 
terms of a Feynman parameter integral. The electric monopole form factor is given by 



4V) = -cm 



M 7 (/x - 7 ) S 



2tt 



1 tan — 

V 4 7, 



g 2 A Mmlj 



(m w + 2 7 ) 



where we have defined the functions 
A(x 



tan 



-l f Q 



\j^i 2 + x(l — x)q 2 /4 



4 7 



3{x) 



i i 

dx — — tan -1 8 
xA 



qx 



(3.12) 



(3.13) 



2( 7 + m 7r + A) ' 

The operator with coefficient D2 does not contribute to these observables. Because of the 
running of C2, the above expression is independent of \i to the order we are working [|TJ. 
From eqs. (|3.10| - |3~T2| ) we determine the charge radius of the deuteron to NLO, 

2\LO 1 



2\NLO 



8 7 z 

MQu- 7 ) 2 + g\Mml{?>m« + 10 7 ) 



167r 7 967r/ 2 7(m 7r + 27) 3 

A comparison with the experimental value is given in §|TV[ 
At NLO, the electric quadrupole form factor is given by 



(3.14) 



dx 



1 



7rl67r/ 2 g 3 Jo x8 4 A 



(3.15) 



x ( [3q V(l + 3 2 ) 2 - 24qm n 3x(l + 3 2 ) + 16m 2 /3 2 (3 + /3 2 )] tan" 1 3 
+3 [-48m 2 /3 2 + 8m w qx3(3 + 23 2 ) - q 2 x 2 (3 + 5/? 2 )]) 
From this expression one can extract the quadrupole moment to first nonvanishing order: 

. nn L " = • — • v..: — • (3.16) 



H'Q 



^ 307r/ 2 (m, + 2 7 )3 

A comparison with the experimental value is given in §|V|. 
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FIG. 3. The expansion ofT 1 , where the photon corresponds to the vector potential A 1 . The cou- 
pling of the photon to the nucleon lines represents the entire one-body current from C\, eq. \2.1\ ), 
including the magnetic moment contribution. The first graph is the LO contribution at 0(Q°), while 
the remaining graphs are the NLO contributions at 0(Q l ). The photon couplings arise through any 
of the operators in Co, C\ or £2- We specifically distinguish the C2, D% and L2 vertices by the 
symbols m, o, and • respectively. 



B. The NLO computation of the magnetic form factor 

In order to calculate the magnetic form factor of the deuteron, we need the matrix element 
of the spatial current (p', k |Jg TO |p, j). This entails computing P, using the coupling of the 
spatial component of the gauge field, Aj, discussed in § [II A| . The expansion of P in Feynman 



graphs is shown to subleading order in Fig (^). Following our power counting rules, T l begins 
at 0(Q°), and so an expansion analogous to eq. (|3.8|) for the matrix element of J* m implies 
that the magnetic form factor has the expansion 



F^ + F^+OiQ 2 ) . (3.17) 



Our task in computing Fm is greatly simplified by recognizing from eq. ( |3.3|) that we need 
only pick out contributions with spin structure anti-symmetric in the deuteron polarization 
vectors. It is straightforward to check that none of the graphs shown in Fig. (Q) contribute 
to Fm when the photon coupling arises from any of the operators N^Y) 2 N, g^N^cr ■ (£D£' — 
£i~D£)N, Tr DfJlD^Yy in eq. ( |2.11| ), or the four-nucleon operator with coefficient Ci in 
eq. ( [2.151) . At LO, only the photon coupling via the isosinglet nucleon magnetic moment 
one-body operator contributes, 



— Korftr ■ BN = ^±^ N^a- ■ BN , (3.18) 
2M u 2 ' v ; 



and we find 



eF ffV) = ± K F (o) ( 2) = 
2M d M c yq ' 




(3.19) 
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For the deuteron magnetic moment this gives /i£P = (/i p +/i n ), simply the sum of the neutron 
and proton magnetic moments. 

At next order, Q 1 , there are contributions to Fm arising from coupling the photon via 
eq. Q3.18 )) along with insertions of the C% operator or one pion exchange; there is also 



a contribution from the two-body current arising from the operator in eq. Q2.15| ) whose 



coefficient is L 2 . We find that there are no pion exchange current contributions at this 
order, nor any two-body current contribution from the C 2 operator in eq. ( |2.15| ). With the 
exception of the two-body contribution involving an explicit factor of B (see eq. ( 2.15|) ), all 



the graphs contributing are all proportional to those giving rise to the electric form factors 
in Fig. (0). Therefore to this order we can express the magnetic form factor in terms of the 
electric form factors and a single new coupling constant. We find 

= Oh + *0 (^V) + V)) + eL 2 l(n - 7 ) 2 , (3.20) 

and the deuteron magnetic moment is given by 

^ i0 = eL 2 ^(/i- 7 ) 2 . (3.21) 

7T 

where L 2 depends on the renormalization scale fi in such a way that is fi- independent. 
A comparison with the experimental value is given in the next section. 



C. Effective range theory 



In effective range theory the electromagnetic form factors are assumed to be dominated 
by the asymptotic S-wave deuteron wave function, 



^(ER)( r ) 



7 



- 7 r 



47r(l — 770) r 



(3.22) 



Assuming the small r part of the deuteron wave function is only important for establishing 
the normalization condition, Fc(0) = 1, the prediction of effective range theory for the form 
factor Fc(q 2 ) follows from the Fourier transform of ^^^(r)! 2 , 



c 



(a 2 



1 



-1 + — tan -1 



7 r o, 



47, 



(3.23) 



This yields the charge radius, 

(r 2 ) ER = 



87 2 1 — 7r 
1 

87^ 



l + 7^o +7 2 ^o + ••• 



(3.24) 



It is instructive to compare the effective range theory prediction for the charge radius with 
that from effective field theory. In effective field theory at NLO the effective range is (tq = 
at LO), 
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r = C 2 (ji) 



M(/i- 7 ) 2 S&M 
2w Anf 2 



8 7 
3 m n 



+ 2 



r 



mi 



(3.25) 



Using this it is straightforward to show that the 7 expansion of the NLO effective field theory 
charge radius, 



(r 2 ) 



57' 
1 



1 + C 2 {n) 



M 7 (/i-7) 2 t g 2 A M-rml(3m w + IO7) 

12^/2(77^ + 2 7 ) 3 



2tt 



+ 



, + C 2 („) + 



4nf 2 



87 7' 



3 



(3.26) 



and the 7 expansion of effective range theory agree to order 7 2 at linear order in tq. 

Effective range theory predicts the matter radius, r m , with remarkable precision. Using 
r = 1.75 fm effective range theory yields r£ R > = 1.98 fm. The most recent measurement 
of the deuteron charge radius is r c h = 2.1303 ± 0.0066 fm from which the matter radius 

(In effective field theory the effects that 



1.9685 ± 0.0049 fm 12 



is found to be r m 

distinguish between the matter and charge radius don't arise until NNLO.) The numerical 
success of the prediction of effective range theory for the matter radius suggests that the 
most important higher order terms in effective field theory are those that arise from iterating 
the NLO potential from C 2 and one-pion exchange. However, from the effective field theory 
perspective this cannot be justified since there are new local operators that will contribute 
at the same order. 

Effective range theory can also be used to predict the magnetic form factor and it gives, 



2M d 



{li p + H n )F { c R \q 



(3.27) 



In the following section Fjf R \q 2 ) and F^ R \q 2 ) are compared with experimental data. 



IV. COMPARISON WITH DATA 

We now compare the analytic results of our effective field theory perturbative expansion 
for the deuteron form factors with experimental data. We have evaluated these expressions 
at the same renormalization point \x = m n used in refs. [|IJ and have used the same value 

C 2 {m n ) = 9.91 fm 4 (4.1) 

derived from a fit to the NN scattering phase shifts in the spin triplet channel. The values 
of Co and D 2 do not enter our expressions explicitly, but they do enter indirectly through 
the constraint on the two-point function that the deuteron pole occurs at the correct binding 
energy, eq. (|A6|) . Given C 2 from the NN phase shift analysis, we have no new parameters 
at through NLO for fitting the electric form factors. As we have seen, for the magnetic form 
factor, a single new parameter, L 2 , enters at NLO. 

We first consider that static moments of the deuteron, at q 2 = 0. We have analytic formu- 
las for the charge radius, the quadrupole moment, and the magnetic moment in eqs. ([3.14 ) 



(ft.16 ), and (|3.21|) respectively. A comparison of these values to experiment is given in 
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TABLE I. Electromagnetic properties of the deuteron 





LO 


NLO 


(LO+NLO) 


Experiment fT^JTH 


RMS charge radius (fm) 


1.53 


0.36 


1.89 


2.1303(66) 


Magnetic moment (N.M.) 


0.88 


-0.02 (fit) 


0.86 (fit) 


0.85741 


Quadrupole moment (fm 2 ) 




0.40 


0.40 


0.2859(3) 



Table Q. The charge radius shows a rapid convergence to the measured value, which is en- 
couraging. The LO calculation is expected to be within ~ 30% of the experimental value, 
while the NLO calculation is expected to be within ~ 10%. It is clear from Table | that 
this expectation is fulfilled. When the NNLO calculation is performed we expect that the 
result is within ~ 3% of the experimental value. The magnetic moment agrees well with 
experiment at LO, and then is fit to the experimental value at NLO by choosing the strength 
L 2 of the two-body magnetic operator appropriately. The LO prediction for the magnetic 
moment is much closer to the experimental value (within ~ 3%) than naively expected from 
the power counting. The quadrupole moment vanishes at LO, and the NLO value of 0.40 fm 2 
is off by ~ 40%, as expected from the power counting. It would be useful to compute the 
NNLO contribution to fiQ to see if it exhibits the same convergence as the charge radius. 
The idea of including pions perturbatively has been used previously to estimate the deuteron 
quadrupole moment |14j], obtaining a value of 0.40 fm 2 . More interesting is that iterated 



potential pion exchange reproduces the deuteron quadrupole moment || reasonably well. 
This suggests that contributions to the quadrupole moment from higher order countert- 
erms are small compared to additional insertions of potential pion exchange. This smallness 
is not something that arises naturally in the effective field theory. It is also interesting 
that state-of-the-art nuclear calculations of the quadrupole moment |15j (~ 0.270 fm ) are 
systematically lower than the experimental value by ~ 7%. This strongly suggests that 
dynamics beyond potential interactions are required, something that effective field theory 
provides a systematic way to include. 

Of greater interest is the comparison of the form factors over a range of q 2 , as we should 
be able to see at what momentum the expansion begins to fail; our naive estimate is that 
the expansion is in powers of q/2A^N ~ g/(600 MeV). The differential cross section for 
elastic electron-deuteron scattering is given by 



da da 
dH = dn 



Mott 



A(q 2 ) + B(q 2 ) tan 2 9/2 , (4.2) 



where A and B are related to the form factors [|TT|j : 

A = F 2 + 2 - V F 2 M + , 

5 = ^7(1 + 7?)^, (4-3) 

with i] = — (p — p') 2 / (4Mj) ~ g 2 /4M|. In order to compare with data, we take our analytic 
results for the form factors and expand the expression eq. ( [4.3| ) in powers of Q, where 
V~0(Q 2 ) 
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100 200 300 400 

q (MeV) 

FIG. 4. A plot of A(q 2 ) vs. q in MeV for elastic electron- deuteron scattering. The dotted curve 
shows the result of the LO calculation, while the solid curve is the NLO prediction. There are no 
free parameters at this order. The dashed curve shows the result of effective range theory. 



A 
B 



24 0) 4 1) 



0{Q 2 



+ 



+ 0(Q 4 



(4.4) 



We see that to the order we are working, A is sensitive only to the electric form factor 
Fc, while B depends only on the magnetic form factor Fm- A comparison of A and B with 
experimental data in Figs. dUU) shows that our expansion is quite successful, and converging 
rapidly, in the kinematic regime where it is expected to work. The data for Fig. [| was taken 
from ref. |l(J, and the error bars are smaller than the size of the points; the data for Fig. |5| 
comes from refs. [ T5 — |T5 1 . 

It is evident from Figs. (HH) that the NLO effective field theory calculation of the 
deuteron form factors in not as accurate as what effective range theory gives. The va- 
lidity of effective range theory over such a wide range of momentum occurs because of the 
smallness of the shape parameter, n. In the effective field theory expansion, the coeffi- 
cients in the effective range expansion themselves have perturbative expansions. However, 
ultimately when carried out to higher orders the effective field theory calculations will be 
more precise than the effective range calculations. This is because the effective field theory 
correctly describes the strong interactions, which effective range theory only approximates. 



V. CONCLUSIONS 

We have demonstrated that one can compute properties of the two nucleon system to 
surprising accuracy simply by calculating several Feynman diagrams. The technique for 
doing this was introduced in refs. where it had been shown how to work at NLO for NN 
phase shifts in both spin singlet and triplet channels. While encouraging, those results were 
not definitive as the NLO calculation required three free parameters in both spin channels. 
The true test of the theory has been presented in this paper with the computation of the 
electromagnetic form factors of the deuteron — by using the parameters fit to scattering 
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0.01 
B(q 2 ) 

0.005 




100 200 300 400 



q (MeV) 

FIG. 5. A plot of B{q 2 ) vs. q in MeV for elastic electron- deuteron scattering. The dotted 
curve shows the result of the LO calculation, while the solid curve is the NLO prediction. There 
is one free parameter at this order, L2 , which is fixed to correctly reproduce the deuteron magnetic 
moment. The dashed curve shows the result of effective range theory. 



data, we are able to reproduce very well at NLO both the electric and magnetic form factors 
in elastic e — d scattering up to momentum transfers q 2 = (400 MeV) 2 = 4.1 fm~ 2 . Since 
our results are analytic, it is straightforward to analyze what features in the data are due to 
short versus long distance physics. A central feature of our expansion — that pion exchange 
is perturbative — is supported by the success of our fit to the form factors. 

One feature of our results which is especially encouraging is the evidence that the expan- 
sion is converging rapidly. This is apparent in the improvement of the fits to e — d scattering 
data in going from LO to NLO, improvements in the static moments of the deuteron. The 
RMS charge radius presented in Table | deviates from the experimental value by ~ 30% 
at leading order, but only ~ 10% at next-to-leading order. The magnetic moment was off 
by ~ 3% at leading order, and exact at next-to-leading order, due to the contribution of 
a new operator. At NLO the results of effective field theory for the electric and magnetic 
form factors of the deuteron, Fc{q 2 ) and i<A/(g 2 ), are not as accurate as those from effective 
range theory. However at NNLO the effective field theory approach should reach (or even 
surpass) the precision of effective range theory. Furthermore, the methods developed in this 
paper can be used to make predictions for other properties of the deuteron, including those 
for which effective range theory is not applicable. 

Since the NLO result for the quadrupole form factor is the first nonvanishing term in 
its expansion, it is expected to work less well. At the level we are working, the quadrupole 
form factor does not contribute to e — d scattering, however, we can compare the quadrupole 
moment with experiment, and it is ~ 40% too large. We expect this error to be substantially 
reduced in the NNLO calculation, which includes among other things the exchange of two 
potential pions, and short distance 3 Si — 3 Di transitions. In general, it would be interesting 
to compare NNLO results for all of the form factors. Other effects that enter at this order 
are relativistic corrections, radiation pions, and nucleon form factors. 

There remain a number of NLO calculations to be done in the two nucleon system, and 
we are optimistic about their success. Extending this procedure to the three body system 



and beyond remains a fascinating challenge [20 
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APPENDIX A: THE GRAPHICAL EXPANSION OF THE MATRIX ELEMENT 

OF J EM 
1. Irreducible Green functions 



In this appendix we derive eq. ( |3.5|) which is central to our calculation of the deuteron 
electromagnetic form factors. We begin with the interpolating field defined in the text, 



N T RN 



(Al) 



where Pi is the projection defined in eq. (|2.14j) . The full propagator G is defined as the time 
ordered product of two of these V fields: 



G(E)5 lJ = / d 4 *e-^-P- x ) (0|T \v\{x)VM] |0> = <% = 



iZ(E) 



E + B + is 



(A2) 



where B is the deuteron binding energy. By Lorentz invariance, the propagator only depends 
on the energy in the center of mass frame, namely 



E = E — 



_P_ 

AM 



+ 



E = (p°-2M), 



(A3) 



where the ellipses refers to relativistic corrections to the dispersion relation. The numerator 
Z in eq. flA"2j) is assumed to be smooth near the deuteron pole, and when evaluated at the 
pole gives the wavefunction renormalization Z, 



Z(-B) 



dG-^E) 
dE 



-i -i 



(A4) 



J E=-B 



It is convenient to define "irreducible" Green functions as the sum of graphs which do not 
fall apart when the graph is cut between incoming and outgoing nucleons at the four-fermion 
vertices proportional to Cq. The irreducible 2-point function is denoted by E, and has the 
expansion shown in Fig. |l|. One can see graphically (Fig. |] that the relation between G and 
E is 



G 



1 + iC S 



It follows that 



E=-B Cn 



1 dE 
l?dE 



i 

Z 



(A5) 



(A6) 
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FIG. 6. The expansion of of the full 2-point function G in terms of the irreducible 2-point 
function E. 




FIG. 7. The expansion of full three-point function G^ in terms of the irreducible two- and 
three-point functions E, T^. 



In general, unphysical quantities such as Z, Co, the deuteron wavef unction, etc. will depend 
on the renormalization scale fi, while S-matrix elements will be /i-independent. 

In order to compute the matrix element of the electromagnetic current between two 
deuteron states, we first define the 3-point function 

G^(E,E',q) = [ d 4 xd 4 ye^ Ex0 -^e^ E 'y°-^ <0|T [vj(x) ^(0)2^)1 |0> , (A7) 



where g M = (E' — E, p' — p) is the photon momentum. G M is related to the desired form 
factor via the LSZ formula 



(p',j \J^\P,i) = Z \G- 1 (E)G- 1 (E')G^(E,E', q 



E, E ^-B 



(A8) 



where G(E) is defined in eq. (|A2|) . It is convenient to reexpress this formula in terms E and 
the irreducible 3-point function, which we call T^. It is easy to see graphically (Fig. [7|) that 
the relation between G M and T M is 



G£ (£X q) 



FUE,E,q) 



l + iC Z(E)) (l + iC Z{E'] 

rfj (e, ~e' ,ol)G (e) g (e' ) 



(A9) 



Making use of this relation and eqs. ( |A5| 
of the current in terms of and E: 

(p',j \Jem\P^) = = Z 



E(£)E(£) 

allows us to reexpress the matrix element 



nm(E') 

dE(£)/d£ 



E,E 



(A10) 



E,E 

It is this relation that has a simple perturbative description in terms of Feynman graphs. 
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2. Computing E 



We can now compute E in our perturbative expansion, writing E as 

CO 

E(£) = £E (n) (E) 



(All) 



n=l 



where E( n )(£') ~ 0(Q n ). The leading contribution to E is shown in the first row of Fig. |T|, 
and is 0{Q) according to the rules of the previous section. These graphs are readily evaluated 
using the formula eq. ( |2.16| ), with the result 



M 

'An 



-ME - is 



(A12) 



The subleading contribution is 0(Q 2 ) and one must compute the three graphs shown in the 
second row of Fig. [l|. The result is Jl[ 



£( 2) (£) 



.g 2 A M 2 mi 



327T 2 / 2 



i tan 



1 /2v / ME\ 1 (ml + AME\ 



1 



-In 

2 



+ 1 



2/2 



C 2 ME + D 2 ml) E (1) (E) 



(A13) 



To the order we are working we truncate the expansion in eq. ( |3.6| ) to the nonrelativistic 
result, 



E~E 



AM 



(A14) 



the first relativistic correction enters at NNLO, or 0(Q 3 ). Other NNLO contributions are 
shown in the third row of Fig. [IJ and include the exchange of two potential pions, or one 
radiative pion (see [J] for discussion) as well as several other graphs. 

From eq. ( |3.5|) we see that what is needed is dH/dE evaluated at E — —B. From 
eqs. ( |AT^ - |AT3l) we find 



dE (1) 



dE 



. M 2 

E=-B ~~ 87T7 



dE 



(2) 



M 3 



E=-B ^16^7 



dE 

where we have defined 



v 2 



ml 



+ 27 



+ D 2 ml (7 - n) - C 2 j{n 27) 



(A15) 



7 



= Vmb 



(Ai6) 
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3. Computing T c 



The leading contribution to the matrix element of the J® m current between deuteron 
states arises from the three-point function r9_j\, the first graph in Fig. ^], 



M-i) 



2nq 



47 



(A17) 



where q — |q| is the magnitude of the photon 3- momentum, and 7 was defined above in 
eq. flXjgp . 

At subleading order we need to sum the diagrams in the second row of Fig. In each 
case, there is a minimally coupled A photon coupled to the proton propagator, with either 
an insertion of the C 2 or D 2 contact interactions, or a single pion exchange^. We find 



1 (o) 



e5i 



M 3 
16tt 2 



+ 7, 



D 2 U) 2^- ^ tan" 1 



^- j + C 2 {fi){fi - 7) ( - 7 - " — tan 



xqA(x 



tan 



-1 



— — tan 1 










))1 


+ 7 + m^) 





+e (q^- - g 2 <5. 
2 



167r 2 / 2 g 3 7o 



x 



f xqr\ 



1 



x 2 q 2 r 2 



sin 



-(A(x)+7+m^-)r 

xr 3 A(x) 
f xqr 



^3 + Sm^-r + m 2 ? 



(A18) 



where A(x) is defined in eq. ( |3.13| ). 

As discussed in the text, the calculation of the parts of T l which are antisymmetric in the 
deuteron polarizations is completely analogous to the complete calculation of T presented 
here. 



APPENDIX B: NO OFF-SHELL AMBIGUITY 

COMPUTATION 



AN EXPLICIT 



When working with potential models for NN interactions one often faces ambiguities 
about how to continue matrix elements off-shell. In an effective field theory approach, there 
is no such ambiguity. All uncertainties arising in a consistent calculation are due to higher 
order operators neglected at the order one is working [JZT 
the effect of the operator 



To illustrate this, we consider 



O = (N T PiN) 



'2M 



-2M 



(Bl) 



4 One might worry that in fact there are four-nucleon contact interactions involving the combina- 
tion a covariant time derivative Dq, and hence a direct photon coupling to the N^N^NN vertex. 
In fact, such an operator may be eliminated by using the equations of motion. We demonstrate 
this by explicit calculation in appendix [B[ 
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where is the gauge covariant derivative^ , 




Dp = d„ + ieQ em A^ , (B2) 

Qem being the electric charge matrix. The operator O is not Galilean invariant but nonethe- 
less we can in principle consider how it enters the NLO calculation of the deuteron form 
factors via the graphs in Fig. ^. However, to the order we are working, it vanishes by the 
equations of motion, 

/ TY2 \ 

iV(x,t) = 0. (B3) 

One might naively think that the equations of motion imply that the operator O will not 
enter a calculation of NN phase shifts (as the nucleons are on-shell in that process), yet 
that O will affect deuteron matrix elements, since the nucleons are not on-shell in a bound 
state. This would mean that a new constant enters the deuteron calculation which cannot 
be determined vis NN scattering. 

However, this is reasoning is incorrect, and we now show by explicit calculation that 
operator O does indeed vanish when considering deuteron matrix elements. This result is 
consistent with general theorems of field theory that state that off-shell matrix elements 
are arbitrary (they can be changed by making a field redefinition) and that the S'-matrix 
elements never depend on them (even when the matrix element is between bound states). 

As an example, consider the contribution to the deuteron three-point function T° of the 
operator O in the graphs of Fig. || corresponding to the matrix element 

T° = <0|T [vj(E, 0)V j {E' 1 q)A°(q , q)] |0> , (B4) 

where E' = E + q°. The first graph, Fig. ||(a), includes the photon-independent part of O 
and a minimally coupled A photon on a nucleon leg. It is proportional to 

r d D ~ l k d°- l \ E-k 2 /M 

(a) K " ! J ( 2.)-. (2.)- {E _ 5) {E : |) {E , : - ■ (B5) 

where D — > 4 at the end of the calculation. To evaluate this integral, we used the fact 
that the first term in the numerator cancels the first propagator, and that in dimensional 
regularization 

d D ~ l k , x 

. (B6) 



(2tt) d - 1 

The second graph in Fig. |8] is similar, and proportional to 



5 To be chirally invariant, the covariant derivative should include pion fields, but as the pion 
couplings do not enter to the order we are working, we have set them to zero. 
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(a) (b) (c) 

FIG. 8. Feynman diagrams contributing to the matrix element of 0\ and Oi- The gray circle 
denotes an insertion of the operator O in eq. (Bl). 



. , cP^k d D - x \ E' - (k 2 + (k + q) 2 ) /2M 

[bj oc 



(27T)"-! (2tt)^ (e - g) - ^g±3l!) - ^IgSl! 

rf^k rf 13 ' 1 ! 1 

(2k)°-i (2tt)^ ( E _ g) ( E , _ lf±lgHf 



(B7) 



Finally, the third graph, Fig. H(c), arises from the A photon coupling in O, and gives 

. r d D - l \z d D ~H 1 

(Cj ^ +Z ■/ W" 1 (27r)^-i (js _ g) {e> - l!±g±al!) ' { ] 

It follows that the sum of the three graphs in Fig. |8| vanishes, and there is no off-shell 
ambiguity arising from this new operator O. Similar remarks hold for other operators 
with a single time derivative. Therefore we can choose to only include the two spatial 
derivative operators in the Lagrangian (our interaction proportional to C2) and to eliminate 
the analogous operators with a time derivative by the equations of motion — even though we 
are considering nucleons bound in a deuteron. This result is not peculiar to the particular 
regularization scheme we used. 



23 



REFERENCES 



[1] D.B. Kaplan, M.J. Savage and M.B. Wise, Phys. Letts. B424 (1998) 390; 

|nucl-th/9802075| , to appear in Nucl. Phys. B. 
[2] A.V. Manohar, Lectures given at 35th Internationale Universitatswochen fuer Kern- unci 

Teilchenphysik, Perturbative and Nonperturbative Aspects of Quantum Field Theory, 

Schladming, Austria, 2-9 Mar 1996. |hep-ph/960622^ . 
[3] C. Ordonez, L. Ray, U. van Kolck, Phys. Rev. C53 (1996) 2086. 
[4] T.-S. Park, K. Kubodera, D.-P. Min, M. Rho, |hep-ph/97ll463 . 
[5] H.A. Bethe, Phys. Rev. 76 (1949) 38. 
[6] H.A. Bethe and C. Longmire, Phys. Rev. 77 (1950) 647. 
[7] T. Ericson, Nucl. Phys. A416 (1984) 281c. 

[8] T. Ericson and M. Rosa-Clot, Ann. Rev. Nucl. Part. Sci. 35 (1985) 271. 

[9] E. Jenkins, M. Luke, A.V. Manohar and M.J. Savage, Phys. Letts. B302 (1993) 482; 

B388 (1996) 866(E), and references therein. 
[10] S. Weinberg , Phys. Lett. B251 (1990) 288; Nucl. Phys. B363 (1991) 3; Phys. Lett. 

B295 (1992) 114. 

[11] See, for example, M. J. Zuilhof, J. A. Tjon, Phys. Rev. C22 (1980) 2369. 
[12] A.J. Buchmann, H. Henning and P.U. Sauer, Few-Body Systems, 21 (1996) 149. 
[13] T. Ericson, W. Weise, Pions and Nuclei, Oxford University Press (1988). 
[14] D.Y. Wong, Phys. Rev. Lett. 2 (1959) 406. 

[15] B.S. Pudliner, V.R. Pandharipande, J. Carlson, S.C. Pieper and R.B. Wiringa, 
|nucl-th/9705009 . 

[16] G. G. Simon, CH. Schmitt, V. H. Walther, Nucl. Phys. A364 (1981) 285. 

[17] B. Grosstete, D. Drickey, P. Lehmann, Phys. Rev. 141 (1966) 1425. 

[18] D. Benaksas, D. Drickey, D. Frerejacque, Phys. Rev. 148 (1966) 1327. 

[19] D. Ganichot, B. Grosstete, D. B. Isabelle, Nucl. Phys. A178 (1972) 545. 

[20] P.F. Bedaque and U. van Kolck, [nucl-th/9710073| ; P.F. Bedaque, H.-W. Hammer and 

U. van Kolck, |nucl-th/9802057| . 
[21] H. D. Politzer, Nucl. Phys. B172 (1980) 349; C. Arzt, Phys. Lett. B342 (1995) 189. 



24 



